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Thermodynamics of itinerant magnets is studied using a classical model with one parameter char- 
acterizing the degree of itinerancy. Monte Carlo simulations for bcc and fee lattices are compared 
with the mean-field approximation and with the Onsager cavity field approximation extended to 
itinerant systems. The qualitative features of thermodynamics are similar to the known results of the 
functional integral method. It is found that magnetic short-range order is weak and almost indepen- 
dent on the degree of itinerancy, and the mean-field approximation describes the thermodynamics 
reasonably well. Ambiguity of the phase space measure for classical models is emphasized. The 
Onsager cavity field method is extended to itinerant systems, which involves the renormalization 
of both the Weiss field and the on-site exchange interaction. The predictions of this approximation 
are in excellent agreement with Monte Carlo results. 



I. INTRODUCTION 

The thermodynamics of magnetic materials is often de- 
scribed using the Heisenberg model in which the spins are 
attached to lattice sites. Real magnets are much more 
complicated, because the magnetization is due to band 
electrons whose degree of localization varies between dif- 
ferent materials. This so-called itinerancy manifests it- 
self in the fluctuation of the magnitudes of the local mo- 
ments, which may be defined in a mufhn tin sphere or us- 
ing a projection in an appropriate basis. Thus, the degree 
of itinerancy may be characterized by the relative impor- 
tance of longitudinal and transverse (rotational) fluctua- 
tions of the local moments.^ In the localized (Heisenberg) 
limit the longitudinal spin fluctuations (LSF) have a large 
energy scale and are suppressed. This limit is approached 
in some magnetic insulators. Metals, on the other hand, 
are often quite far from this limit, because the exchange 
splitting and the bandwidth are typically of the same or- 
der. Experimentally, itinerancy is most clearly revealed 
in the paramagnetic susceptibility by the deviation of the 
effective moment found from the Curie- Weiss constant 
from the true local moment, as well as by the deviations 
from the Curie- Weiss law. 

A large amount of work has been devoted to the ther- 
modynamics of itinerant magnets using phenomenologi- 
cal Ginzburg-Landau models for weak ferroniagnetsi'^'^ 
or the Hubbard model and the functional integral 
methodsiii^i^i^ These studies have clarified the role of 
LSF in thermodynamics and explained the observed be- 
havior of the paramagnetic susceptibility. However, these 
methods are unsuitable for quantitative studies of re- 
alistic materials. Ginzburg-Landau expansions, as is 
well known, correctly describe only the contribution of 
long-wave fluctuations and must always be rigged with a 
wavevector cut-off. Such models are useful in the stud- 
ies of critical phenomena, but they are irrelevant to the 
determination of the critical temperature itself, which 
is determined by short-range fluctuations.^ An unsatis- 
factory signature of Ginzburg-Landau models is the ab- 
sence of any information on the short-wave components 



of the exchange interaction in the resulting expressions 
for the Curie temperaturei^>^!^ In our opinion, the ne- 
glect of short-wave fluctuations in these models makes 
their predictions for magnetic short-range order (MSRO) 
also unreliable. The functional integral method, on the 
other hand, suffers from the necessity to make severe and 
ambiguous approximations.- 

Magnetic thermodynamics has also been studied using 
density functional theory (DFT) by treating spin fluc- 
tuations within the adiabatic approximation— assum- 
ing that the relevant fluctuations are well represented 
by constrainedii noncoUinear ground states. The most 
widespread approach is the disordered local moment 
(DLM) approximatior^iSii^ which relies on the single-site 
approximation and is designed to approximate the DFT 
ground state of a system with random directions of the 
local moments. The LSF have been neglected in all im- 
plementations of this approach so far, restricting its ap- 
plication to magnets which are close to the localized limit. 
In particular, the DLM method neglecting LSF fails for 
(strongly itinerant) nickel where it finds vanishing local 
moment in the paramagnetic phase. 

Other authors studied itinerant thermodynamics by 
mapping the results of first-principles energies for var- 
ious spin configurations (including both transverse and 
longitudinal fluctuations) to a classical Hamiltonian in 
which variable local moments play the role of dynam- 
ical variables, and then exploring the thermodynamics 
of this Hamiltonian using either the variational principle 
in reciprocal spacei^ or Monte Carlo (MC) simulations 
in real spacc i^^i^^'^^ These calculations clearly show that 
LSF, as expected, are very important in nickel. More- 
over, they revealed only weak magnetic short-range or- 
der (MSRO) above the Curie temperature for both Fe 
and Ni, which is similar to the Heisenberg model. These 
results are consistent with the fact that in any lattice 
model with no frustration all correlation corrections to 
the mean-field approximation (outside of the critical re- 
gion) should be small in the parameter where z is the 
number of neighbors within the interaction rangeJ^ On 
the other hand, very strong MSRO above Tc was foundi^ 
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in Ni using the ab initio spin dynamics method, which, 
similar to DLM, is based on the adiabatic approximation 
and neglects LSF. 

Classical models with variable local moments seem to 
capture the important qualitative features of the thermo- 
dynamics of itinerant magnets which are similar to the 
predictions of the functional integral method. However, 
these models have been built and studied only for a few 
particular materials, and a general study of their ther- 
modynamic properties has not been undertaken. Such a 
study is useful as a step to more refined models with the 
advantage that numerically exact results for a classical 
model are easily accessible through Monte Carlo simula- 
tions. Therefore, in this paper we explore the thermody- 
namics of a classical spin fluctuation model as a function 
of the degree of itinerancy using MC simulations and sim- 
ple analytic approximations. We emphasize that here we 
are not concerned with the "mapping" procedure (which 
can be quite challenging) but rather focus on the other 
separate part of the program, i.e. on the determination of 
thermodynamics once the Hamiltonian has been defined. 
We therefore restrict ourselves to the simplest possible 
realization of this model which includes only one free 
parameter characterizing the degree of itinerancy. 



II. MODEL 

Our model is a lattice version of the phenomenological 
model of Murata and Doniacb^ written with a vector 
order parameter—: 

H= i^x"^(q)mqm_q-|- j^mf 
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Here rrii denotes the magnetic moment at site i whose 
length is unrestricted, and / the Stoner exchange- 
correlation parameter. We have separately written the 



local 



dBi/dmi and nonlocal J^. 
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of the unenhanced inverse susceptibility. This model in- 
volves a number of simplifying assumptions: (1) It is 
classical in the sense that nij are dynamical variables and 
not operators. (2) Both local and nonlocal parts of the 
inverse susceptibility are considered to be independent 
of the magnetic state and isotropic. In general, xTj^ is 
a Cartesian tensor which depends on the magnetic state 
and reduces to a scalar only in the paramagnetic state. 
(3) Nonlinear effects are included only through a local 
fourth-order term, similar to the Murata-Doniach model. 

Model ll]) is somewhat similar to that used to repre- 
sent the unified spin fluctuation theory^ classically (see 
Ref. 11!, Ch. 7, and also Ref. [20h . with an important 
difference: the energy of LSFs is included as a function 
of local dynamical variables rrii, rather than that of one 
global parameter (mf). This difference is similar to that 



between the Heisenberg model and the spherical approx- 
imation to it. 

In the ground state all local moments are parallel and 
we recover the Stoner model which is ferromagnetic if 
IN{Ef) > 1, where N{Ef) = x(0) is the density of 
states at the Fermi level in the nonmagnetic state. This 
Stoner criterion can also be written as (/ -f Jq) > Xqq^ 
where Jo — -^v- On the other hand, in the paramag- 
netic or non-magnetic matrix, local moments exist in the 
Anderson sense only if / > Xoo which is stricter than the 
Stoner criterion. We will call this the Anderson criterion. 
(Note that Xoo ^ Vxoo-) 

Introducing reduced local moments ~ nii/mo, 
where mo is the value of all at T = 0, the Hamil- 
tonian ([T]) can be conveniently parameterized: 
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where E{x) — [ax'^/2 + bx^/4:]/Jo with a = Xqo ^ ^ and 
b = Buiq = Jo — a. For the nearest neighbor model with 
coordination number z we have Jnn/Jo = 1/^, and for 
the given lattice H' contains only one parameter, which 
we define as a = arctan5/a. Note that & > is equivalent 
to the Stoner criterion, and a < is equivalent to the 
Anderson criterioni^i 

To understand the meaning of the parameter a, con- 
sider the ground state of Hamiltonian H with a single- 
site excitation, i.e. the state with m; = mo for all i 
except i = c. The energy of this state has a minimum at 
nif, = mo and its curvature with respect to the longitudi- 
nal fluctuation of m^ is = Jo + 26, while the curvature 
with respect to transverse fluctuations is K± — Jo ■ Their 
ratio K^^/K± = 1 + {2b/ Jo) characterizes the relative im- 
portance of longitudinal and transverse fluctuations. If 
6 3> Jo, the fluctuations are mainly transverse, and we 
have the localized (Heisenberg) limit for which a k, ~b 
and a « 37r/4. If 6 ^ Jo, the transverse and longitudinal 
spin fluctuations are equally important; this limit corre- 
sponds to a = 0. The Anderson criterion is equivalent to 
a > 7r/2. Thus, the parameter a characterizes the degree 
of itinerancy and is similar to those appearing in other 
theoriesii^. Note that we always have K\JK±^ > 1, even 
though the macroscopic longitudinal stiffness is propor- 
tional to b and tends to zero at a ^ 0. 

Evaluation of the thermodynamic properties involves 
taking a trace over the quantum states, or a functional 
integral over the classical degrees of freedom. To our 
knowledge, in all classical models reported so far and 
based on ab initio calculations, the uniform measure in 
the space of m^ was usedJ^ii^ii^iii However, our dynam- 
ical variables are not canonical, and therefore the phase 
space measure (PSM) is not known. In the case when 
LSF are important, the PSM has to be supphed along 
with the Hamiltonian as an additional phenomenologi- 
cal ingredient. Strictly speaking, it is not possible to 
disentangle the measure from quantum statistics; for ex- 
ample, in the atomic limit only integer moments with 
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atomic multiplet degeneracies should be present. Ambi- 
guity of PSM is intrinsic to all microscopic classical spin 
fluctuation models including the classical version of the 
"unified theory" of Moriya and Takahashi (Ref. ff], Sec. 
7) and its extensions, as well as the functional inte- 
gral approach combined with the static approximation 
which destroys the correct quantum operator properties. 
In the latter case, the Hubbard-Stratonovich transforma- 
tion can be applied with the interaction term written in 
different ways, which produce different results after the 
static approximation is madei^i^ Two particular choices 
discussed by Hubbard^ result in different measures in the 
space of fluctuating fields Vji uniform in one case, and 
involving the weighting factor v^"^ in another. To ex- 
plore the influence of PSM on thermodynamics, we will 
consider these two measures in the space of the local mo- 
ments va.i. 

III. THERMODYNAMIC PROPERTIES: 
MONTE CARLO AND MEAN-FIELD RESULTS 

Monte Carlo simulations for model ^ were performed 
using the Metropolis algorithm for bcc and fee lattices 
with nearest neighbor exchange. At each step the new 
random direction and magnitude of the moment on one 
site was tried, and sampling of the moment magnitude 
was performed according to the chosen PSM. We used 
supercells with up to 3456 or 6912 sites for bcc or fee lat- 
tices (12 X 12 X 12 unit cells with periodic boundary con- 
ditions). The reduced Curie temperature tc = Tc/ {Jotuq) 
was found using the fourth-order cumulant method^^, 
and the paramagnetic susceptibility was calculated us- 
ing the fluctuation-dissipation theorem. 

In the mean-field approximation (MFA) the magneti- 
zation is found from the self-consistency condition (xz) = 
d\nZi/d{Phw), where 

Z, = r (3) 
Jo 

is the single-site partition sum, hw = (xz) the reduced 
Weiss field, and g{x) the weighting factor, which is either 
1 or for the two chosen PSM's. E{x) is defined after 
Eq. ([2), and f3 = 1/t is the inverse reduced temperature. 

Fig. [T] shows the temperature dependence of magneti- 
zation, the average square of the local moment and the 
paramagnetic susceptibility using the reduced variables 
according to Eq. Results are shown for two val- 

ues of a: 0.487r and 0.697r. In both cases the agreement 
between MC and MFA results is very good for all proper- 
ties (MFA overestimates Tc by 20% or less). The results 
strongly depend on PSM, especially in the more itiner- 
ant case a = 0.487r. In particular, for the uniform PSM 
a second-order phase transition occurs for both values of 
a, but for the PSM with g{x) — x^"^ the phase transition 
is of first order for a — 0.487r, and Tc is nearly 2.8 times 
smaller compared to that for g{x) = 1. 
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FIG. 1; (Color online) (a-b) Reduced magnetization {x^), (c- 
d) mean squared local moment {x^), and (e-f) inverse param- 
agnetic susceptibility a function of the reduced tem- 
perature t = T/{Jorn^). MFA results are shown by solid 
(blue online) lines for g{x) = 1 and by dashed black lines for 
g{x) — x~^. MC results are displayed by black circles for 
g{x) — 1 and by red (gray) squares for g{x) = x~^ (in both 
cases the symbols are filled for fee and empty for bcc lattice). 
The inset in panel (e) highlights the region close to tc for the 
bcc lattice with g{x) = 1 and also shows the results of the 
generalized Onsager method (black line connecting the MC 
points). 



As seen in Fig. [TJ below Tc the average (a;^) declines 
with temperature due to the decrease of the Weiss field, 
which causes the maximum of the distribution function 
to shift to smaller moments. This is in agreement with 
earlier resultsii^i^ii^'i^iii. The width of the distribu- 
tion function increases with temperature, which coun- 
teracts the decrease of the local moment. The PSM with 
g{x) = x~^ puts less weight on the states with large mo- 
ments, and hence (a;^) drops much faster compared to the 
uniform PSM. If the Anderson criterion is not satisfied 
(a < 7r/2) then the most probable moment in the para- 
magnetic state is zero. In this case, (a;^) increases with 
temperature above Tc as seen in Fig. [TJ;. On the other 
hand, if the Anderson criterion is satisfied, the local mo- 
ment may slightly decrease in a range of temperatures 
above Tc, as seen for g{x) — x^'^ in Fig. [TJi. 

The magnetic susceptibility above Tc is shown in Figs. 
le,f. In MC simulations it is calculated using fiuctuation- 
dissipation theorem, while in MFA we directly consider 
the response of the system to the external magnetic field. 
Excellent agreement between MFA and MC is observed 
except for the small error in Tc- In MFA one obtains 
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(4) 



This formula looks similar to the Curie- Weiss expression 
in the Heisenberg model, but here (x^) depends on tem- 
perature, which leads to a renormalization of the CW 
constant and deviations from the CW law. The CW con- 
stant C = dx~^ /dt (for a second-order phase transition) 
is now given by 



1 - 



dlog(x^) 
dlogt 
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Thus, in addition to the usual Heisenberg term the Curie 
constant has a contribution due to the temperature de- 
pendence of (x^) (second term in square brackets in ([5])). 
As a result, the effective moment squared x'^g = 3/C 
deviates from (x^). As discussed above, (x^) usually in- 
creases with temperature above Tc, which, according to 
Eq. ([3]), reduces C and increases xl^. Moreover, for 
the uniform PSM (x^) increases faster with temperature 
compared to PSM with g{x) = x~^, and hence the CW 
constant is much smaller in this case (see Fig. If and also 
le, where the transition is however of first order). 

In Fig. [2] some thermodynamic properties of the sys- 
tem are plotted as a function of the itinerancy parameter 
a. From Eq. ^ it follows that the MFA value of tc for 
the second-order phase transision is found by solving the 
equation 3tc — (x^(tc)), where (x^(t)) is fully determined 
by E(x) in Eq. This is an easy way to estimate Tc for 
an itinerant system using first-principles data for E{x), 
JorriQ, and the assumed PSM. However, for PSM with 
g{x) = x~^ the transition is of first order except for a 
small region close to the local moment limit (in MFA the 
tricritical point where the order of the phase transition 
changes is at atr = 0.6327r). Therefore, in general one 
must consider the minima of the free energy as a func- 
tion of the magnetization, which can also be easily done 
in MFA. Note that the order of the phase transition de- 
pends on the details of the model and can change if, for 
example, the dependence of the exchange parameter on 
the magnetization is taken into account. In particular, 
the phase transition for the model of Ni is of first order 
in Ref. [ij (as seen from the abrupt drop of M{T) and 
Ms at Tc in their Fig. 2) and in Ref. [Ij (as seen from 
the abrupt drop of m in their Fig. 6), even though the 
uniform PSM was used in both of these models. 

From Fig. [5] we see that when the transition is of sec- 
ond order, MFA overestimates Tc by about 20%, which 
is typical for the Heisenberg model. When the transi- 
tion is of first order, MFA gives an almost exact Tc- It 
is important that even for the second-order transition 
the overestimation of Tc in MFA does not depend on 
the degree of itinerancy. This is consistent with the fact 
that the degree of MSRO, which is shown in Fig. for 
T = l.lTc, is quite small and stays essentially constant in 



FIG. 2: (Color online) (a) Reduced Curie temperature tc and 
(b) MSRO parameter {cosd„n) ai T = l.lTc as a function 
of the itinerancy parameter a for the bcc lattice. Solid black 
line, red (gray) squares, and blue (dark gray) circles show the 
results of MFA, MC, and the generalized Onsager method for 
g{x) = 1, respectively. Dashed black line and empty black 
squares depict MFA and MC results for g{x) = x~^. Green 
(light gray) triangles represent the incomplete Onsager reac- 
tion field correction with the on-site interaction left unrenor- 
malized. The blue (gray) dash-dotted line in the upper panel 
shows the effective moment a;^jf found from the Curie con- 
stant for g{x) = 1 in MFA. Very similar results were obtained 
for the fee lattice (not shown). 



the whole range of a. Thus, in our model itinerancy does 
not lead to strong short-range order. This result agrees 
with Refs. ITsIITtI where weak short-range order was found 
for the models of Fe and Ni. Note that if the exchange 
interaction extends to more than one shell of neighbors 
and stays mainly ferromagnetic, the MFA validity cri- 
terion is satisfied even better, and the MSRO parameter 
should further decrease. Similar to the Heisenberg model, 
strong MSRO may only be expected in low-coordinated 
lattices or in the presence of frustration when for some 
pairs Jij/kTc is not small. 

The square of the effective moment x^g is also shown 
in Fig. [5] for the uniform PSM (dash-dotted line). In the 
local limit Xcff naturally tends to 1. However, as a is de- 
creased towards zero, the ratio x^g/ (x^(tc)) increases and 
eventually becomes much larger than 1. Similar behavior 
is found in functional integral theories^ 



IV. GENERALIZED ONSAGER CORRECTION 
FOR ITINERANT SYSTEMS 

Onsager introduced the concept of a cavity field in the 
theory of polar liquids, which is designed to go beyond 
the molecular field approximation (MFA) by including 
short-range order effectsi^ The cavity field is the effec- 
tive internal field which orients polar molecules in the 
ferroelectric phase. Onsager observed that each molecule 
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polarizes the surrounding liquid and thereby generates a 
reaction field acting back on the molecule. However, this 
field is always parallel to the molecule's dipole moment 
and hence does not affect its orientation. Therefore, for 
a liquid with permanent dipoles the reaction field must 
be subtracted from the mean molecular field, the result 
being the cavity field. Onsager also noted that the reac- 
tion field enhances the dipole moments of real molecules 
due to their polarizability. 

The cavity field method was successfully applied to 
Ising2i and Heisenberg2S, magnets which have perma- 
nent magnetic moments. Cyrot^^ noted that Moriya- 
Kawabata's self-consistent renormalization theory for the 
Hubbard model may be essentially reproduced by using 
Onsager-like arguments; more recently this method was 
implemented numerically.-- However, the actual physics 
there is very different; Cyrot's approach seeks the corre- 
lation correction with respect to the Hartree-Fock solu- 
tion, which is unrelated to short-range order. Onsager's 
method was also applied to itinerant nickel^i^, but, as 
we will see below, correct generalization to itinerant sys- 
tems with LSF re quir es an additional ingredient which 
was missed in Ref. [l3l. 

We now generalize Onsager's method to magnets with 
LSF described by Hamiltonian ([1]). Consider model ([T]) 
above Tc in a small external coUinear magnetic field 
Hf^^Bz- We pick site and integrate out the degrees 
of freedom from all the other sites in the partition func- 
tion to obtain the effective Hamiltonian in the form of a 
generating functional for the lattice with a cavity Ex- 
panding this functional around the atomic limit to order 
1/^ we obtain 

= E{m,) mo + E ^o.(m^)^ 



(6) 



where the superscript c refers to the lattice with a 
cavity, i.e. with site removed, and we used the 
fiuctuation-dissipation theorem to express the pair cor- 
relator through the susceptibility. 

In order to find the magnetization and susceptibil- 
ity of the lattice with a cavity we need to solve the 
"impurity problem." Using the linked-cluster expansion 
technique)^ the longitudinal susceptibility of the original 
lattice can be written as follows: 



X = H + HW^H 



(7) 



where W^^s the effective interaction that satisfies the 
equation W = J + JUW, and H is the 1-bond-irreducible 
"polarization operator" which may be shown to be local 
to first order in (All quantities in Eq. ([7]) are ma- 

trices in site indices.) Removal of site may be formally 
represented by a perturbation AH = —TlooSoiSoj to H. 
(The renormalization of Hjj for j ^ due to removal 



of site is at least of order l/z^.) Thus, denoting the 
effective interaction matrix for the cavity lattice as Wc, 
we may write - W^^ = -AH. Using ^ and the 

fact that H is diagonal, we find 



Xi 



Xi3 



Xoo 



(8) 



The average local moments = (m^) for the lattice 
with a cavity are: 
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ext ^ - 2^ Mo, 



Xoo 



(9) 



where Mi are the average local moments of the complete 
lattice without the cavity. The value of Hg^* does not 
affect (as expected), therefore in the right-hand side 
of ^ we may take Mi and Mo for the actual field distri- 
bution. 

From the effective Hamiltonian ^ we can find the 
magnetization at site 0: 



where 



w 



Mo - x°Hw 



Xm_ 

Xoo 



Mo 



(10) 



(11) 



is the renormalized effective field (cavity field), and x° is 
the renormalized bare (atomic-limit) susceptibility. The 
latter may be written as — ('7i^)a/3T, where the av- 
erage paramagnetic squared local moment {m'^)\ is cal- 
culated using a renormalized on-site exchange / = / -I- A 
with A = '^ij JoiJojXij- This renormalization of the bare 
susceptibility is the essential ingredient needed to extend 
Onsager's theory to itinerant magnets. It has no effect 
in the localized limit where is constant. 

As usual, we now obtain the Fourier transform of the 
susceptibility: 



1 - x°( Jq - A) ' 



(12) 



where A — X]q -^qXq/Xoo- We used the same symbol A 

as above in the definition of /, because these expressions 
are identical, as can now be shown with the help of Eqs. 
Cll) and dll). Eq. ^ with the definitions of A, x° and 
/ form a closed set of equations for the paramagnetic 
susceptibility. Note that (jl2p automatically leads to a 
sum rule xoo = x''; which agrees with the fluctuation- 
dissipation theorem. 

At the Curie temperature Xq diverges at q = 0. There- 
fore, from we obtain Tc = ^Jo{m'^{Tc))x/G, where 
G = X)q(l ~ Jq/Jo)~^ is the diagonal element of the lat- 
tice Green's functioni^ Note that the value of A at Tc is 
equal to Jo(l — G"-^) and independent of the degree of 
itinerancy a. 
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The reduced Curie temperature tc and MSRO param- 
eter (cos^nn) at T = l.lTj, calculated in this way are 
shown in Fig. [2] for the bcc lattice and the PSM with 
g{x) — 1. The agreement with MC results is excellent in 
the whole range of a. We repeated these calculations for 
the fee lattice and found excellent agreement with MC as 
well. The accuracy of the predicted tc may be seen from 
Table 1. Similar performance for bcc and fee lattices sug- 
gests that this approximation is not very sensitive to the 
connectivity of the lattice. The paramagnetic suscepti- 
bility is also shown in Fig. [T^ for a = 0.487r, bcc lattice, 
and uniform PSM. The agreement with MC results is 
essentially perfect outside of the narrow critical region. 

TABLE I: Reduced Curie temperature tc for bcc and fee 
lattices for PSM with g{x) — 1: Results of the mean-field 
approximation (MFA), Horwitz-Callen approximation (HC), 
generalized Onsager method (GO) and Monte Carlo (MC). 



a/TV 


bcc 


fee 




MFA 


HC 


GO 


MC 


MFA 


HC 


GO 


MC 


0.032 


0.621 


0.449 


0.451 


0.462(1) 


0.621 


0.465 


0.466 


0.480(2) 


0.148 


0.660 


0.484 


0.486 


0.504(2) 


0.660 


0.501 


0.502 


0.520(5) 


0.250 


0.681 


0.503 


0.504 


0.525(2) 


0.681 


0.519 


0.520 


0.540(5) 


0.352 


0.699 


0.518 


0.520 


0.543(2) 


0.699 


0.535 


0.536 


0.562(2) 


0.422 


0.712 


0.529 


0.530 


0.553(1) 


0.712 


0.546 


0.547 


0.570(5) 


0.483 


0.723 


0.539 


0.541 


0.568(1) 


0.723 


0.557 


0.558 


0.584(2) 


0.553 


0.745 


0.555 


0.557 


0.585(2) 


0.745 


0.572 


0.574 


0.600(1) 


0.602 


0.765 


0.570 


0.573 


0.600(2) 


0.765 


0.589 


0.590 


0.617(2) 


0.687 


0.834 


0.619 


0.622 


0.654(3) 


0.834 


0.640 


0.642 


0.672(6) 


0.735 


0.942 


0.683 


0.688 


0.732(2) 


0.942 


0.708 


0.711 


0.753(6) 


0.750 


1 


0.713 


0.718 


0.770^^ 


1 


0.740 


0.743 


0.788(3) 



The first-order terms in the l/z expansion derived 
above introduce two corrections to MFA. The first one 
is the subtracted mean reaction field; this correction re- 
duces the magnetization. This is the only correction in 
Onsager's method for systems with permanent moments. 
The second correction described by the last term in Eq. 
([6|) adds back the fluctuating reaction field which is al- 
ways parallel to the moment on the central site. For the 
Heisenberg (or Ising) model this second correction has 
no effect, but in itinerant systems it always increases the 
local moments and hence the Curie temperature. There 
is a strong cancelation between these two corrections in 
itinerant systems, and improvement compared to MFA 
may be achieved only if both of them are included. In- 
deed, if the renormalization of the Stoner parameter is 
not taken into account (i.e. if the last term in Eq. ^ is 
dropped) , we find a spurious strong suppression of Tc for 
itinerant systems, as shown in Fig. [2^. 

It is interesting to compare the generalized On- 
sager method with the Horwitz-Callen (HC) approxima- 
tion which is based on the "ring subset" of diagrams 
for the generating functional $ in the linked-cluster 
technique. In this method, the second-order self-field 
G2 is found by differentiating <I> with respect to the renor- 



malized second cumulant M2 , while M2 is represented by 
an integral containing G2 as a parameter. This technique 
does not assume any particular form for the atomic limit, 
and therefore it can be used in our case including LSF as 
well. In the HC method, the on-site correlator may be 
found as i^oo — M2 -A- 2M|G'2, and the sum rule Xoo = 1 
is not satisfied in the paramagnetic Heisenberg magnet. 
However, it is easy to check that the value of Kqq at Tc 
is smaller than 1 by less than a percent in bcc and fee 
lattices. In Onsager's method for the Heisenberg model, 
the sum rule Kqq = 1 is used to fix M2 instead of the in- 
tegral representation as in the HC method. The results 
for Tc are therefore very close. We found that this close 
similarity remains in the entire range of a, as seen from 
Table 1. The generalized Onsager's method is, however, 
technically much simpler. 



V. CONCLUSIONS 

We have studied the thermodynamics of a simple clas- 
sical spin fluctuation model allowing for a variable degree 
of itinerancy. This model is qualitatively similar to those 
used before to study the thermodynamics of Fe and Ni 
using first-principles data j^^'^^i^^ It is worth emphasizing 
that the main drawback of using classical spin models of 
this type is the ambiguity of the phase space measure. As 
we showed above, the thermodynamics is very sensitive to 
this measure for systems with even intermediate degree 
of itinerancy. While the energetics of constrained spin 
configurations may, at least in principle, be accurately 
mapped using DFT calculations, it is not known (to our 
knowledge) how and whether the phase space measure 
can be supplied in a realistic way. 

In the present work, we focused on the general fea- 
tures of the model rather than on the determination 
of its parameters from principles. We found that the 
thermodynamic properties are similar to the results of 
the functional integral approach/i^^i^ Further, we found 
that the mean-field approximation is qualitatively valid, 
and short-range order is weak and almost independent 
on the degree of itinerancy up to the strongly itinerant 
limit where the paramagnetic susceptibility is dominated 
by longitudinal fluctuations. This is in agreement with 
earlier results for the models of Fe and Ni;^^'^^ it is clear 
that this is a general feature of the classical model with 
no frustration. 

Further, we generalized the Onsager cavity field 
method to itinerant systems using an expansion around 
the atomic limit to first order in l/z. Both the inter- 
atomic exchange constant and the Stoner parameter are 
renormalized by short-range order. When both these cor- 
rections are included, the Curie temperature is in excel- 
lent agreement with Monte Carlo results. However, sim- 
ple subtraction of the Onsager reaction field is a very 
poor approximation. 
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